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Abstract. Particle motion is considered in incompressible 
two-dimensional flows consisting of a steady background 
gyre on which an unsteady wave-like perturbation is super- 
imposed. A dynamical systems point of view that exploits 
the action-angle formalism is adopted. It is argued and 
demonstrated numerically that for a large class of problems 
one expects to observe a mixed phase space, i.e., the occur- 
rence of "regular islands" in an otherwise "chaotic sea." This 
leads to patchiness in the evolution of passive tracer distri- 
butions. Also, it is argued and demonstrated numerically 
that particle trajectory stability is largely controlled by the 
background flow: trajectory instability, quantified by various 
measures of the "degree of chaos," increases on average with 
increasing |do/d/|, where (o(I) is the angular frequency of 
the trajectory in the background flow and / is the action. 



1 Introduction 

This paper deals with the kinematics of fluid particles in un- 
steady incompressible flows on the Cartesian plane. Namely, 
we study properties of trajectories (x(t) 7 y(t)) that satisfy 
equations of the form 



i = d y y/, y = -d x y, 



(la) 



where the overdot stands for time derivative and \j/(x,y,t) is 
the streamfunction. Furthermore, we consider the latter to be 
split into a steady background component and an unsteady 
perturbation component, i.e., 



\lf=\lf i0 ^x,y) + ew il) (x,y,t), 



(lb) 



where e is a dimensionless parameter. Equations (1) consti- 
tute a canonical Hamiltonian system with \j/ the Hamiltonian 
and (x,y) the generalized coordinate-conjugate momentum 
pair. 
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Two related issues are addressed in this paper. First, 
we investigate a cause of "patchiness" in passive tracer 
distributions, i.e., distributions that are mostly vigorously 
stirred but include poorly stirred regions (Sect. 3). Second, 
we study the influence of the background flow on particle 
trajectory stability (Sect. 4). Prior to discussing these issues, 
the kinematic models that we use to illustrate our results are 
briefly described in Sect. 2. The conclusions of the paper are 
given in Sect. 5. 



2 Kinematic models 

Two background flow structures in a region [0,L] x [0,W] of 
the j3 plane are considered here. One is chosen to represent 
a large-scale single-gyre wind-driven ocean circulation with 
streamfunction given by (Stommel, 1966) 
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e b - L )/(e b + L -e b - L ), and 



where a := tW/(kXD), b : = (1 

2. Here, D is the 
depth, T the wind stress amplitude (per unit density), and 
X the bottom friction. The other background streamfunction 
chosen corresponds to solid body rotation, 



R: y 



(0) 



COr 

2 



(x-L/2) 2 + (y-W/2) 2 



The reason for this highly idealized choice will be discussed 
below. Parameter values used in our numerical work are 
listed in Table 1 . 

The perturbation streamfunction is constructed by su- 
perposing standing Rossby-like modes with a power-law 
spectrum, namely, 

V/' 1 ' =a^Ae~ r *sin(foc+0jsin(/>'+</) / )cos(c7f + </> fJ ), (2) 
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Fig. 1. Background flow streamlines (dashed lines) along with streamlines corresponding to a snapshot of the total flow at t ss 9 y (solid lines). Panels a 
and b correspond to background flow S with different wave-like perturbation fields superimposed; panel c corresponds to background flow R with the same 
perturbation that was used to produce panel a. 



Table 1. Background flow parameters. 



Parameter 


Value 


L 


10 Mm 


W 


27tMm 


D 


200 m 


8 


9.8 ms- 2 


./() 


10~ 4 s~' 


P 


10"" m^V 


z 


2 x 10~ 3 m~V 2 


X 


10~ 5 s" 1 


(Or 


27ly- [ 



where 

Bk 

a{k,l):= ~e 

and the <j>(k,l)'$ are random numbers uniformly distributed 
between and 2n. Here, Lk/n and Wl/n are positive inte- 
gers; y is a constant; /q is the reference Coriolis parameter; 
and g is the acceleration of gravity. 

Dashed lines in Figs. la,bandFig. lc are streamlines for 
background flows S and R, respectively. Solid lines in these 
figures are total flow streamlines corresponding to a snapshot 
of the flow at t » 9 y. The perturbation in each case involves 
10x10= 100 modes. In Figs. la,c the perturbation has e = 
0.05 and y = 0. In Fig. lb the amplitude of the Lk/n = 1 = 
Wl/n mode is set to zero, (j) k = = <j> { so the flow vanishes 
at the boundary, e = 0.25, and y = 0.4 Mm -1 . 

The flows used to produce Fig. 1 and all of the numer- 
ical particle trajectory simulations presented in this paper 
were chosen to illustrate important aspects of Lagrangian 
dynamics; the flows are in many ways not representative of 
realistic oceanic flows. We note, however, that we focus on 
flows with complicated time dependence, and that the strong 
perturbations to the background are considered. In Fig. lb, 
for example, it is seen that the perturbation leads to the 
presence of an eddy-like structure in the flow. Also, we note 
that in the flows that we have described, particle trajectories 
are periodic in the limit of zero perturbation strength with 



typical periods of about 1 y. Thus in an integration time 
of 10 y most trajectories will have made approximately 
10 revolutions around the gyre. The phenomena described 
below are not limited to gyre-scale flows. In general, 
the trends that we describe should be evident after times 
in excess of a few periods of particle revolution in any 
background gyre flow on which a perturbation field with a 
broad band of frequencies is superimposed. 

3 Passive tracer patchiness 

In this section we present numerical evidence and a theoret- 
ical argument that suggest that for a large class of systems 
of the form (1) phase space (x,y) should be partitioned into 
"regular islands" in a "chaotic sea." Such a mixed phase 
space leads to patchiness in passive tracer distributions. 
Numerical results are presented for a time-periodic flow 
[n = 1 term in the sum in (2)] and subsequently for flows 
with complicated time dependence (n large). 

Figure 2 shows, for the time-periodic case (n = 1), a 
Poincare section and, in the same environment, two ad- 
ditional trajectory diagnostics whose applicability is not 
restricted to time-periodic flows. The Poincare section 
was constructed by plotting the (x,y) coordinates of several 
trajectories at integer multiples of the period of the stream- 
function; it shows the usual mixture of "regular islands" in an 
otherwise "chaotic sea" (cf., e.g., Tabor, 1989). The middle 
panel shows, for a dense set of trajectories with x(0) = xq 
fixed and y(0) = yo variable, a plot of y vs. yo at a fixed value 
of t. The initial conditions chosen fall inside the region of 
the Poincare section shown, and it is seen that both regular 
islands and the chaotic sea evident in the Poincare section 
can be identified in the y vs. yo plot. The same structures 
can also be seen in the lower panel of Fig. 2 which shows, 
for the same trajectories used to produce the middle panel, 
finite time estimates of Lyapunov exponents (described in 
more detail below), v vs. yo. Plots of y vs. yo and v vs. yo 
are used below to distinguish between apparently regular and 
apparently chaotic trajectories for flows with complicated (« 
large) time dependence i.e., in flows for which a Poincare 
section cannot be constructed. 
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Fig. 2. Poincare section (top), final vs. initial meridional position (middle), 
and finite-time estimate of the Lyapunov exponent as a function of initial 
meridional position (bottom). The middle and lower plots were constructed 
by tracking 10 4 particles (Ayo ss 5.5 km) for a duration of 150 y in 
background flow S with a time-periodic perturbation superimposed. Particle 
initial positions fall on the horizontal line shown in the top panel. The 
perturbation parameters chosen were e = 0.015 ,7 = 0, kL/n = 3 = IW/x, 
<j>'s = 0, and 27t/a = 0.25 y. 



Figure 3 shows plots of y vs. yo and v vs. yo for the 
nonperiodic flows used to produce Fig. 1 . Trajectories in Fig. 
3b are generally more unstable than in Fig. 3a. The enhanced 
stability in Fig. 3a is reflected in a relatively unstructured 
y(yo) plot and smaller (on average) Lyapunov exponents than 
are seen in Fig. 3b. In both cases the background flow 
structure is the same; the difference in the stability behavior 
is due to the difference in the strength of the perturbation. 
As expected, trajectory instability is seen to increase with 
increasing perturbation strength. 

The difference in trajectory stability seen in Figs. 3a,c has 
a different explanation. The same perturbation was used in 
both cases, so this cannot be the cause of the difference. The 
cause is the influence of the background flow; this topic will 
be discussed in detail in the following section. 

We return our attention now to Fig. 3b which corresponds 
to the strongly perturbed flow shown in Fig. lb. It is 
seen in Fig. 3b that embedded among mostly chaotic 
trajectories are bands of apparently nonchaotic trajectories. 
These nonchaotic bands are most readily identified among 
the trajectories whose initial positions are near the center of 
the gyre; the reason for this will be discussed in the following 
section. Bands of nonchaotic trajectories far from the gyre 
center are also present, however. This is seen in Fig. 4 
where two regions of Fig. 3b are blown up. These apparently 
nonchaotic bands of trajectories are the counterparts of the 
"regular islands" seen in Fig. 2. Trajectories in these 



bands diverge only very slowly (power law dependence on 
time) from neighboring trajectories while chaotic trajectories 
diverge at an exponential rate from neighboring trajectories. 
The nonchaotic regions of flows are important in applications 
because they correspond to regions where the concentration 
of a passive tracer will remain high for a long duration. The 
existence of these regions leads to a large variance of tracer 
concentration or "tracer patchiness" (cf., e.g., Pasmanter, 
1988; Malhotra et al., 1998). 

Another way to visualize passive tracer patchiness is 
offered in Figs. 5 and 6. Both figures show discrete samples 
of a material line of fluid at t = (vertical line segments in 
the figures) and at t = 12 y, in the environments shown in Fig. 
1. Initial conditions in Fig. 5 are as in Fig. 3, whereas those 
in Fig. 6 are as in Fig. 4. These figures again show that while 
most of the initial material line segment is vigorously stirred, 
there are small portions of the initial segment, corresponding 
to the island-like structures seen in Figs. 3 and 4, that are 
poorly stirred. 

We turn our attention now to explaining the occurrence of 
island-like structures in Figs. 2-4. First, we note that in the 
background flow, particle motion is describable using action- 
angle variables, reviewed below, and trajectories fall on tori. 
For perturbed systems with periodic time dependence, as in 
Fig. 2, it is well-known that particle trajectory dynamics 
are constrained by the KAM theorem (cf., e.g., Arnold, 
1989) which guarantees that for sufficiently small e some 
of the original tori — and associated nonchaotic motion — 
are preserved. Related theoretical results, generally known 
as KAM theory (cf., e.g., Tabor, 1989), describe how the 
nonsurviving tori break up to form chains of "islands" 
surrounded by a "chaotic sea" as seen in Fig. 2. For a 
large perturbation strength £ all of the original tori will have 
been broken up, but the secondary islands that are formed in 
the process are robust and persist even when the magnitude 
of the perturbation exceeds that of the background flow. It 
has been shown (Brown, 1998; Beigie et al., 1991) that for 
multiply-periodic perturbations the situation is essentially 
the same as for perturbations with simple periodic time 
dependence. This follows from the observation that (1), with 
y/{x,y, ait,--- ,<7„f) where a it is defined modulo 2n, can 
be transformed to an autonomous Hamiltonian system with a 
bounded phase space with (n + 1) degrees of freedom that is 
constrained by n integrals. KAM theory (the KAM theorem 
and related results) applies to the transformed system, so 
phase space is generically partitioned into nonintersecting 
regular and chaotic regions. A Poincare section could, 
in principle, be constructed for such a system by using a 
multiple slicing technique (cf. Parker and Chua, 1989) but 
slicing is practical only when n = 1. The significance of 
the extension of KAM theory to multiply-periodic systems 
is that in the system defined by (1) with y/ 1 ' given by (2), 
phase space (x,y) is expected to be partitioned into "regular 
islands" in an otherwise "chaotic sea." The numerical 
evidence presented in Figs. 3-6 supports this expectation. 

The coexistence of regular and chaotic fluid particle tra- 
jectories in mesoscale and large-scale oceanic flows has been 
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Fig. 3. Trajectory final meridional position y (top panels), finite-time estimate of the Lyapunov exponent v (middle panels), and absolute value of stability 
parameter a (bottom panel) as a function of initial meridional position. Background and perturbation fields are as in Fig. 1. All particles have initial 
longitudinal positions xq at the center of the background gyre. The integration time is 12 y and Ayg ~ 5.5 km. 
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Fig. 4. Blow up of two portions of the top and middle panels of Fig. 3b using Ayo ~5m. 



suggested in some analyses of surface drifters and submerged 
floats (Osborne et al., 1986, 1989; Richardson et al., 1989; 
Brown and Smith, 1990). The preceding discussion provides 
an explanation of the underlying physics. 

4 Particle trajectory stability 

In this section we describe the important influence of the 
background flow on particle trajectory stability that was 
mentioned above in our discussion of Figs. 3a,c. The ideas 
presented here apply to any canonical Hamiltonian system 
in which the Hamiltonian consists of a superposition of 
an integrable component and a nonintegrable perturbation. 
Other applications are described in Beron-Vera and Brown 
(2003a,b). 

The explanation of this behavior makes use of the action- 



angle description of the motion of particles in the background 
flow (cf., e.g., Abdullaev and Zaslavsky, 1991). Let 

/:= ^^dxF(x;v/ (0) ), 

where Y is the meridional coordinate of an isoline of 
be the action variable, and consider the canonical transfor- 
mation (x,y) i— > (#,/), defined implicitly by 

y = d x G, # = <? 7 G, G(x,I):= j dxY(x; y/ (0) ), 

where & is the angle variable. According to the above 
transformation, 

W{x,y,t)^ v/(°)(/)+£V/«(/,iV) 

and Eqs. (la) take the form 

/ = -£<?,> y/ (1) , ?J = co + e<3/y/ (1) , (3) 
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Fig. 5. Initial position (straight segment) and position after 1 2 y of a discretely sampled (Ayo ~ 
plus perturbation) flows of Fig. 1. Colors indicate the initial meridional position of particles. 
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Fig. 6. Blow up of two portions of Fig. 5b corresponding to trajectories shown in Figs. 4a,b; here Ayo ~ 5 m. Colors indicate the initial meridional position 
of particles; the initial material lines are so short that they appear as dots. 



where 



a(I) :=dy/(°Vd/ 



(4) 



When e = 0, Eqs. (3), which have one degree of freedom, 
are autonomous and the corresponding Hamiltonian, 
is an integral of motion that constrains the dynamics. As a 
consequence, the equations can be solved by quadratures and 
the motion is periodic with angular frequency CO. Namely, 
I = Iq and # = #o + cot mod27T, where Iq and #o are 
constants. Every solution curve is thus a line that winds 
around an invariant one-dimensional torus {Iq} x T 1 cRx 
T , whose representation in (x,y)-space is the closed curve 
given by the isoline y/ ' = \jf^°\lo). 

With the perturbation term, the corresponding Hamilto- 



nian, yrw + exjf ' , is no longer an integral of motion (the 
equations are nonautonomous) and the system may be sensi- 
tive to initial conditions, thereby leading to chaotic motion. 
The distinction between regular and chaotic trajectories is 
commonly quantified by the Lyapunov exponent (cf., e.g., 
Parker and Chua, 1989), a measure of the rate at which 
neighboring trajectories diverge, 



Voo := lim — In I v Q | 



(5) 



where v Q (f) is the largest of the two eigenvalues of the so- 



called stability matrix Q(f), which is given by 



Q 



dk 1 
d k d d 





' o" 








(**)" 


<o' 


{si) 


+£ 





Because of the area preservation property of Eqs. (la) or 
(3) the product of the two eigenvalues of Q is unity, so 
there is no loss of generality in considering only the largest 
eigenvalue. Each column of Q corresponds to a vector 
perturbation (5/, 5t?) to a trajectory in the nonautonomous 
system (3), and satisfies the so-called variational equations, 

81 

(6) 

where co' := dco/dl. Equations (6) and (3) constitute a 
system of four coupled equations. 

Variational equations that describe the growth of pertur- 
bations using Cartesian coordinates (8x, 8y) have the same 
form as (6) except that in the Cartesian form all four elements 
of the first matrix on the r.h.s. of (6) are generally nonzero. 
Our numerical finite-time Lyapunov exponent estimates are 
based on the Cartesian equivalent of Eqs. (6) and (3), which 
is generally more convenient for numerical calculations. We 
have chosen to show the (81, 5$) form of these equations to 
highlight the important role played by co' . An example of 
a closely related study which does not exploit action-angle 
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Fig. 7. Evolution of a material line of fluid in background flows S (a) and R (b) 
a function of initial meridional position. 

variables, and which consequently overlooks the critical 
importance of to', is Richards et al. (1995). 

A simple but very important observation follows from 
the action-angle formalism. Dependence of both particle 
(3) and variational (6) equations on the background flow 
enters only through the function co(I). Equations (6) strongly 
suggests that trajectory stability and to' are closely linked. 
The following heuristic argument explains the mechanism by 
which to' is expected to control trajectory stability. If one 
assumes that £ is small and the second derivatives of yrW 
are zero-mean random variables, then when ft)' = these 
terms should lead to slow (power-law) growth of 8$ and 
81. If [ft)' | is large, this term will cause a rapid growth of 
5$ | for any nonzero |5/|. The perturbation terms will then 
lead to a mixing of \8~&\ and \Sl\. The term ft)' will lead, 
in turn, to further growth of \8~&\. As this process repeats 
itself, both |5/| and |5#| are expected to grow rapidly. 
The role played by ft)' in this process is to amplify small 
perturbations caused by the second term on the r.h.s. of 
Eqs. (6). Thus when e is small, trajectory instability is 
expected to be significantly enhanced when |ft)'| is large. 
When e is sufficiently large that the two terms on the r.h.s. 
of Eqs. (6) have comparable magnitude, the role played 
by | ft)' | in amplifying perturbations is expected to be much 
less important. Increased trajectory instability should result 
in larger numerical estimates of Lyapunov exponents. A 
dynamical-systems-based argument on the role of to' in 
controlling trajectory stability is given below; that argument 
is consistent with the above heuristic argument. 

The lower panels of Fig. 3 show the absolute value of the 
stability parameter (Zaslavsky, 1998; Beron-Vera and Brown, 
2003a,b) 

as a function of trajectory initial condition; recall that these 
initial conditions correspond to variable yo with xq fixed at 
the gyre center. Comparison of the middle and lower panels 
of Fig. 3a suggests that when the perturbation to the back- 
ground steady flow is weak, trajectory instability increases, 
on average, with increasing \a\. Figure 3b shows that for 
a strong perturbation this trend is less strong, although the 
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most stable trajectories are clearly those in the region of 
the flow where \a\ is small. The background flow R used 
to produce Fig. 3c was chosen because it has the property 
a = V/. Because the same perturbation flows were used to 
produce Fig. 3a and 3c the difference between these figures 
is entirely due to the difference in the background flows. The 
remarkable stability of trajectories in Fig. 3c is due to the 
property a = V/ in flow R. The same comment applies to 
the difference between Fig. 5a and 5c, which were produced 
using the same flows that were used to produce Fig. 3a and 
3c. All of the aforementioned observations relating to Figs. 
3 and 5 are consistent with the heuristic argument given in 
the preceding paragraph describing how |ft)'| is expected to 
control trajectory stability. 

The physical interpretation of the stability parameter a 
is illustrated in Fig. 7. Figures 7a,b show the evolution 
of a material line of fluid in background flows S and R, 
respectively. The material line is shown at t = and at t = 12 
y. Also shown is a plot of both 2 n/ to and \a\ as afunctionof 
yo (for yo > yc, the meridional coordinate of the center gyre) 
in each environment. As a consequence of the uniqueness 
of solutions to Eq. (1) and continuity of the velocity field, 
the material line of fluid cannot break or intersect itself but 
it can increase in complexity with time. Because the motion 
in Fig. 7 is integrable (i.e., each point of the material line is 
constrained to lie on a surface of constant / Vf) and because 
attention is restricted to background flows for which has 
compact and closed level sets, i.e., gyre flows, the length of 
the material line can grow with time, at most, following a 
power law. Background flow R has a special property. In that 
background flow the material line just rotates clockwise at a 
constant rate to = Or (= 2% y _1 ), independent of /, so a = 
V/. In contrast, to varies with / in background flow S. The 
monotonic decay of to as a function of / in background flow 
S induces a shear in phase space which causes the outermost 
points of the material line to rotate more slowly than the 
innermost ones and, hence, causes the material line to spiral. 
In background flow R there is no shear. In polar coordinates 
radial shear can be defined as 

rd r {r- l u e ), (8) 
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where r is the radial coordinate and uq is the 0-component 
of the velocity field. More correctly, this quantity is twice 
the rf?-component of the strain-rate tensor for rotational 
motion (cf., e.g., Batchelor, 1964). The connection with 
motion in phase space can be accomplished by identifying 
/ with r and col with ug. The replacements r \— ► I and 
uq i ► col in (8) thus give the analogous expression Ico' for 
the shear in phase space. Notice that this expression is 
(apart from the o -1 -factor) the stability parameter a. We 
have chosen to include the co~ l -factor in the definition of 
a because of precedent (Zaslavsky, 1998; Beron-Vera and 
Brown, 2003 a,b) and because it is convenient to make a 
dimensionless. 

To see the importance of the shear in the background 
flow, compare Figs. 7a,b with Figs. 5a,c, which show the 
evolution of the same initial material line segments in the 
total (background plus perturbation) flows. Notice the highly 
complicated structure of the segment in the perturbed flow S 
(Fig. 5a) as compared to that in the unperturbed one (Fig. 
7a). (Note that the number of particles used to produce Fig. 
5a is far too small to resolve what should be an unbroken 
smooth curve which does not intersect itself.) In contrast, 
observe that in environment R the perturbation has only a 
very minor effect on the evolution of the material line (Figs. 
5 c and 7b). 

Additional insight into why a should be expected to 
control trajectory stability comes from the following argu- 
ment. The perturbation streamfunction ei// 1 ' has the effect 
of introducing perturbations to the action / of a given particle 
by the amount 81. If 81 is assumed to be small and of the 
same order as the perturbation streamfunction [Si = 0(e), 
say], then CO experiences the change 

co i — > (1 + a8l/I)co 

+ 0(e 2 ). The perturbation to CO depends on both the per- 
turbation 8 1 and the background flow via a. Under the 
change / i— ► / + 5/, a sufficient condition for co to remain 
invariant at 0(e) is a = 0. This provides an explanation 
for the remarkable stability of the particle trajectories in 
flow R. To 0(e) a nonvanishing shear (a ^ 0) appears as a 
necessary condition to sustain the successive stretching and 
folding of the material line of fluid after it gets distorted by 
the perturbation. (Of course, chaotic motion is still possible 
when a — provided that e is sufficiently large.) It is thus 
expected that where \a\ is small (resp., large) there will be 
less (resp., more) sensitivity to initial conditions and, hence, 
the motion be more regular (resp., chaotic). Support for this 
conjecture is given in the numerical simulations presented in 
this paper. 

Finally, the role of co' in dynamical systems theory de- 
serves further comment. A nondegeneracy condition, co' ^ 0, 
must be satisfied in order for the KAM theorem to apply 
and, hence, to guarantee that some trajectories are nonchaotic 
provided the strength of the time-dependent perturbation is 
sufficiently weak. This theorem does not imply, however, 
that trajectories are unstable when CO 1 = 0; the KAM theorem 



does not address this limit. The mechanism that leads to 
chaos is the excitation of resonances at discrete frequen- 
cies. For a sufficiently strong perturbation, neighboring 
resonances overlap and chaotic motion results (cf., e.g., 
Tabor, 1989). The width in frequency of each resonance is 
proportional to Ico'] 1 / 2 , so one expects, on average, motion 
to become increasingly chaotic as \co'\ increases. This ex- 
pected trend is consistent with the arguments and numerical 
simulations that we have presented above. The trend toward 
increasingly chaotic motion with increasing \co'\ does not, 
of course, rule out some nonchaotic motion for fixed but 
large co' | . Note also that the trends that we have described 
apply on average; details depend on details of the flow, both 
background and perturbation. 

5 Concluding remarks 

In this paper we considered particle motion in unsteady 
incompressible two-dimensional flows consisting of a steady 
background gyre on which a highly structured unsteady 
wave-like perturbation is superimposed. The numerical 
simulations presented strongly suggest that: (i) phase space 
is mixed, characteristic of near-integrable one-and-a-half- 
degree-of-freedom Hamiltonian systems; and (ii) particle 
trajectory stability strongly depends on the structure of the 
background (steady) component of the flow. 

The mixed phase space structure, in which "islands" 
of stability emerge from an otherwise chaotic "sea," was 
explained as a consequence of the applicability of KAM 
theory. The mixed phase space provides an explanation 
for the occurrence of patches of poorly stirred fluid in a 
mostly vigorously stirred flow. Trajectory instability was 
shown to increase with increasing magnitude of a := Ico' / co, 
where 2n/co(I) is the period of revolution of a particle in the 
background gyre flow and / is the particle's action variable 
in the background flow. 

These results provide important insight into the physics 
underlying Lagrangian ocean dynamics. In addition to this 
insight, the results described are potentially important in a 
variety of practical problems. The occurrence of Lagrangian 
"islands of stability" has important implications for the 
transport and dispersal of tracers ranging from nutrients to 
toxic pollulants. Knowledge that such "islands" are smaller 
and less abundant, on average, in regions of flows where \co'\ 
is large might be exploited when deciding where to place a 
sewage outfall, for example. 
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